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I— I , Abstract 

u, 

■ It is shown that optimal network plans can be obtained as a limit of point 



allocations. These problems are obtained by minimizing the mass transportation 
on the set of atomic measures of prescribed number of atoms. 



\o 1 Introduction 
(N 

^ ! Optimal mass transportation was introduced by Monge some 200 years ago and is, 

OO 

O ' today, a source of a large number of results in analysis, geometry and convexity. 
' 

O ■ Optimal Network Theory was recently developed. It can be formulated in terms of 



Monge-transport corresponding to some non-standard metrics. For updated references 
on optimal networks via mass transportation see [BS, BCM]. 

In this paper we restrict ourselves to the transport of a finite number of points. 
Consider points {xi, . . . xn} (sources) in a state space (say, M.''), and another 
points {yi, . . . yj\f} C M''' (sinks). For each source Xi we attribute a certain amount of 
mass rrii > 0. Similarly, m* > is the capacity attributed to the sink yi, while 

N N 
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We denote this system by an atomic measure A := — A where 

N N 

\+ = ^mi6^^ ; X- = ^m*6y^ (1.1) 
1 1 

where is the Dirac deha function. 

The object is to transport the masses from the sources to sinks in an optimal 
way, such that the sinks are filled up according to their capacity. A natural cost was 
suggested by Xia [X]: For each g > 1, W^'^\X) is defined below: 

Definition 1.1. Given A as in U.l\) . 

1. An oriented, weighted graph (7, m) associated with X is a graph 7 embedded in M.^ , 
composed of vertices V{'^) and edges E{^i). The orientation of an edge e G £'(7) 
is determined by de = v'^ — v~ where vf G V{'^) are the vertices composing the 
end points of e. The graph 7 and the capacity function m : -£^(7) U {0} 

satisfy 

(a) {xi, . . . xtv, 1/1, . . . vn} C V{j). 

(b) For each i e {1,N}, J2{e,x,ed+e}^e = rrii and J2{e,y,ed-e}^e = m*, where 



± 



(C) For each V G V{j) - {Xi, . . .y^}, E{e;^,Ga+e} "^e = E{e;«G9-e} 

2. The set of all weighted graphs associated with A is denoted by r(A). 
3. 

W^''\X) := inf V lelmy" (1.2) 

{7,m)Gr(A) ^ 

There are two special cases which should be noted. In the limit case g = 1 the op- 
timal graph satisfies V{j) = {xi, ...Vn} and W'^^^X) = Wi{X+, A"). Here Wq{X+, A") 
for g > 1 is the Wasserstein distance between A+ to A~, 



N N 
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Wg{X^, A ) := ( mill ^ XI ~ y^\'^''i 



the minimum is taken in the set of N x N matrices satisfying 

N N 

> , ^ 7i J- = m* ■ ^ 7i,j = rrii . 
i=i j=i 

In particular, Wi depends only on the difference A = — A~ (which is not the case 
for g > 1). 

The second case is the limit q = oo. This is the celebrated Steiner Tree Problem 
[HRW]: 

inf > lei , 

7Gr(A) ^ 
eGE(7) 

where, this time, r(A) is the set of all graphs satisfying {xi,yj;mi,m* > 0} C ^(7) 
and is, actually, independent of the masses rrii and capacities m* (assumed positive). 

In [W, Thm 2] it was shown that Wi is obtained from Wq by an asymptotic expres- 
sion for the limit of infinite mass: 

Theorem 1.1. //A = A+ — A~ is any Borel measure satisfying J dX = 0, then 
lim M^-^/'? min W„ (/i + A+, /i + A") = l^i(A+,A-) 

where B\j stands for the set of all positive Borel measures fi normalized by J dfi = M. 

If, in particular, A is an atomic measure of the form (11.11) . than it can be shown 
that for fixed M the minimizer of Wq [fi + , fi + X~) in is an atomic measure of 
a finite number of atoms as well. 

The main result of the current paper demonstrates that the network cost W'^'^^ is 
obtained by similar expression, where the total mass M is replaced by the cardinality 
of the support of the atomic measure fi. 
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2 Main results 



For each n G N, let be the set of all atomic, positive measures of at most n atoms, 
that is: 

I j=l 

Theorem 2.1. For any q> 1 and A as in / fi.ij) 

lim n^-i/'? inf W^, (/i + A+, + A") = I^^^VA) . (2.1) 
The set S"*"'" is, evidently, not a compact one. Still we claim 
Lemma 2.1. For each n G N, a minimizer yU„ G 

W^"\a+,A-):= inf l^J/i + A+,/i + A-) =l^g(/i„ + A+,/x„ + A-) (2.2) 

Remark 2.1. Note that W^^'' depends on each of the component A^ while the limit 
p^('j) = lim.n^oon^^^^'^W^^^ depends only on the difference A = A+ — A^. 



Theorem 2.2. Let fin be a regulafQ minimizer of Wq {jj, + A"*", /i + A~) in S"*"'". Then 
the associated optimal plan spans a reduced weighted irecl {'^ni^n) which converges (in 
Hausdorff metric) to an optimal graph (7, m) G r(A) of ^.2) as n ^ 00, 

3 Auxiliary results 

We first reformulate W , as given by (12. 2p . in terms of a linear programming: 



^see Definition 13.11 
See Definitions 13.21 and 
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Given q > 1, n eN, Z = {zn+i, . . . zn+u) e A = A+ - A" as given by ([H]) 

and 7 := {7^^ 1 <ij <n + N} e T{n, A+, A") : = 

n+N 



7ij > , 1 < A; < iV ^ ^ 7fc, 



n+N 



i=l 

n+7V n+N 



i=l 



Yl = Yl ^^'^ iV + 1 < J < n + iV I , (3.1) 



i=l 



Let 



where 



n+N n+N 



F,(Z,7):=5^^7m^m(^) 



1 1 



We observe 



\Zi - 


1 


1 




- 








+,A- 





1 < i,j < N 



inf F,(Z,7) 

ZG(R'=)",7Gr(n,A+,A-) 



(3.2) 



(3.3) 



(3.4) 



Our first object is to prove Lemma 12.11 that is, to replace the "inf" in (13. 4p by 



"min". 



Definition 3.1. 7 G r(n, A^,A ) is called a regular plan if it satisfies the following 
for any 1 < i, j < n + N: 



(a) if k > 1 and ii = i,i2, ■ ■ - ik = i then LI^^i IhM+i ~ 0- T-^'"- po-rticular 7jj7j,j = 

and -fi^i = 0). 

(b) If k > 1, k' > 1 and {ii = z, Z2, • • • 4 = j} ^ {'^'1 = 4; • • • 4' ~ 

If 'J is a regular plan, then fi G B^'"' is called a regular measure if for each i G 
{A^ + 1, . . . n + A^} there exists z-i G M.^ where ^{{zi}) = Yl^=i 



Lemma 3.1. For each Z G (M'^)" and any plan 7 G r(n, A"*", A ) there exists a regular 
plan Y ^ r(n, A^, A^) satisfying Fq{Z,Y', ) < Fq{'Z,l)- 

Proof, a) Assume Ti^iZllii^ii^^ > 0. Let ii^ such that 7i,(,,i,o+i < 7i,,i,+i for any 1 < / < 
/c. Then Tl'^'i^+i := Ik^k+i - 7j,(,,i,o+i while Yij = lij otherwise. It follows that 
e r(n, A+, A") and Fg{Z,Y') < Fg{Z,Y- Thus Y' verifies (a). 

b) We may assume that {^2, • • • ik-i} ^ ih^ ■ ■ ■ \' _i} = ^ for, otherwise, choose 2 pairs 
of indices ii = i[, and im = i'^' for which {ii+i, . . . im-i} n . . . J = 0. 

Assume (^'1=111^,11+1^ ^n|LY^77., ^ > 0. Assume (with no limitation to gen- 
erality) that Yl^^il^k - > YA=i\^i[ ~ ^i!,,!^^^- Let ii^ such that 
^ any 1 < / < fc. Then set 

Y' ' ■■= 7? / +7''' • 
^^+1 ^A+i ^''0'%- 



while = 7[^- otherwise. Then Y verifies (13. ip while 



F,(Z,7l=F,(Z,7-)-7[;,,„+i 



fc-l k -1 

<F,(Z,7''^)<F,(Z,7) (3.5) 

□ 

Lemma 3.2. The set of regular plans in B^''^ associated with r(ra, A"'",A~) ^3.2^ is 
compact. 

Proof. Let Zi be some point in the support of yU where ix{{zi}) = Q. We show an a- 
priori bound on Q (hence compactness). By (13. ip there exists a point where •ji^i^ > 
Q/{N + n). We can define such a chain i = ii,i2, . . . where 7j;,i;+i > fl{{zi^}) / {n + N) . 
In particular it follows that fi{{zi^}) > Q/{n + NY^^. By part (a) of the definition of 



regular plans, this chain must be of length ar most n. By (13.1 p it must end at some 
i;- := j G {1, . . . A^}. So, li{{zj}) > Q/{n + N)''-'^ > Q/{n + NY'^ . On the other 
hand, < maxi<K7v max{mi, m^*} := M so Q < M{n + A^)"""^. □ 

Corollary 3.1. For fixed Z E (R^)", A satisfying ( fi.i)) and q > I, the function Fg 
admits a minimizer 7 G r(n. A"*", A~). Moreover, this minimizer is regular. 

Proof, of lemma 12.11 

For a fixed A satisfying fll.ip and g > 1 it follows from Corollary 13.11 that 

F,(Z,7):= inf Fg{Z,j)= min F,(Z,7). 

7er(n,A+,A-) ^ 7Gr(n,A+,A-) ^ 

It is also evident that Fg is continuous and coercive on (M'^)" and that W^^\x~^, A") = 
min2g(iRfc)n Fg{Z, A). In particular (13. 4p is attained at a pair (2', 7) where 

n+iV 

/^n := ^^'^^^^ ^ ^""'^ (3.6) 

is a regular minimizer of (12. 2p . □ 

Next we associate a weighted graph (7,m) with a transport plan 7 G r(n, A"^, A~) 
and Z G (M^')" as follows (see Fig 1) 

Definition 3.2. Let Z = {-Zat+i, . . . z^^n} ^ (M'^)" anc? 7 G r(n, A+, A^) . The associ- 
ated weighed, directed graph (7, m) is defined as 

(1) = {xu...yN,Zi,...Zn} := {Cu ■ ■ ■ Cn+2N} ■ 

(a) E{^) is given by the set of segments ek,i ■= [Ck,Ci] for which 'jk,i > 0, while 
dek,i = Cfc - 0- 

(Hi) m^^ ^ := 7fcj. 

(iv) For each Zi G ^(7), deg{zi) := 7^^ + 7^,^ > 0}. 
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Figure 1: A tree associated with a regular transport plan (A^ = 2, n = 11 



Lemma 3.3. Let {Z,'~f) as in Definition \3.2\ where i is a regular plan inr{X,n). Then 
the associated graph (7,m) contains no cycles. In addition, \E{'y)\ < n + 2N^. 

Proof. The result that the graph 7 contains no cycles follows directly from Defini- 
tion 13.11 -a. 

It follows that any vertex v G V^(7) must belong to a chain Cij := Ci, • • • Cfc where 
k < n, (i = Xi and (k = yj. By Definition 13. li b there exists at most one such a chain 
for any pair (xj, yj) G {xi, . . .xn} x {yi, . . . yN}- In particular there exists at most N'^ 
such chains. 

Let now Q G Cij. If the degree of (i is greater than 2, there exist deg{Q) — 1 > 1 
chains which contain Q. By Definition 13 .li b it follows that if two chains Ci' j', 
intersect the chain Cij then either C-' = Cf (and, in particular, they intersect Cij 
at the same point), or i" 7^ i' and j" 7^ /. Hence the number of chains crossing Cij is 
bounded by 2N. As the number of chains {Cij} is bounded by A^^ it follows that there 
exists at most 2A^^ chains which intersect other chains. Hence '^y^y(^Xf-j{deg{v) — 2) < 
2A^^ which implies the result. □ 

Next, we elaborate some properties of an optimal regular plan. 

Definition 3.3. A chain of a regular plan is a sequence of indices ii, . . . ,ik such that 
> Z^'" k > I > 1 while 'j^j = for any j G {1 . . . , n + 2A^}. A maximal chain 
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is a chain which is not contained in a larger chain. 

Remark 3.1. By (Ejp we also get that 'yii^i,,^,^ is a constant along any maximal chain 
ii, . . . ,ik where 1 < I < k. 

Lemma 3.4. If'y is a regular optimaA plan then for any chain {Cii, • • • Ofc}? Ci;+i ~Cii = 
(i^, — (i^, for any l,l' G {1, ... A; — 1}. In particular, all points on a chain of the 
associated directed graph corresponding to an optimal plan are equally spaced on a line 
segment . 

Proof. If is a regular optimal plan then Z = {zi, . . . Zn) is a minimizer of Fq(Z, •Jq) 
in (M'^)". In particular ^ = holds for any I < j < n. If j = ii is embedded in a 
chain then by definition and Remark 13.11 we obtain 



1^. — ^. 1 9 2 — 7- \9 2 







which implies the result. □ 

Let us now re-define the associated directed graph (7, m) corresponding to an op- 
timal regular plan (see Fig 2) 

Definition 3.4. The reduced weighted graph (7^, m) associated with an optimal regular 
plan is obtained from (7, m) (Definition \3.2^) by identifying all edges corresponding to 
a maximal chain {ii, . . .i^} with a single edge [ChXik] ^'^^ assigning the the common 
weight = to this edge (recall Remark WH]) . 

Corollary 3.2. A reduced weighted graph (7^,m) associated with an optimal regular 
plane in B^''"' satisfies the following: 

(i) All the vertices of are of degree at least 3. 

(ii) The number of vertices of is at most 2N^ where N is the number of atoms of 

(in particular, independent of n). 



Figure 2: The reduced version of the tree presented in Fig 1: All vertices of degree 2 
removed. 

(in) All the edges of are line segments. 

(iv) There exists C > 0, depending only on N, such that C > rrif. > 1/C for any 

(v) There is a compact set K C which contains 7^ for any n G N. 

Proof. Part (i) follows directly from Definition I3.4[ Part (ii) from Lemma 13.31 part 
(iii) from Lemma 13.41 To prove part (iv) we repeat the proof of Lemma 13.21 with the 
additional information of (ii) (that is, the bound on the number of edges is independent 
of n). Part (v) is evident. □ 



4 Proof of Theorems 2.1 and 2.2 



Proof, of theorem 12.11 

Let (7, m) be a weighted graph. Then by the Holder inequality 



1/9 



eg£;(7) 



^ m\'^\e\ < me|e|^ 1^(7) l^""^)/"? 



(4.1) 
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If, moreover, (7, m) is obtained from a regular plan 7 G r(n, A+, A ) then 

iy^«(/i„ + A+,/i„ + A-)< ^e|er (4.2) 

where /i„ G S"^'" associated with 7 via fl3.6p . By Lemma [2.11 there exists an optimal 
measure fin £ iS"^'". Hence (14.21) holds with an equality for this choice of Moreover, 
fin can be chosen to be a regular measure (Definition 13.11) hence, by fl4.1|4.2p and by 
Lemma 13.31 

Wq{X)< J2 ^y'|e| < W^"\A+,A-)|n + 2Ar3|('?-i)/^ . 

eeE{y) 

This implies the inequality 

liminf n^-^/wi"^ (A+.A") > W^^^X) . 

To prove the reverse inequality in (12.11) we consider an optimal weighed graph (7,m) 
of Wq{X) and construct yU„ G i?"*"'" supported on 7 which satisfy 

lim n^-i/Wg + A+,/i„ + A-) = V m^'^lel = 1?('?)(A) . 

ee£;(7) 

Assume ne is the number of points of /i„ on the edge e, and any atom of fin in e is of 
weight rrie. The contribution to W^ji^fin + X^,fin + A^) from e is, then 

|e|\'' me\e\'^ 

Up 



^ 0-1 

riej nl 
n^-'Wq^{fin + A+, fin + X-) ^ n'^-i J] 

The constraint on rig is given by ^eg^;^^) ?T-e = 't-- Let us rescale We '■= ne/n. Then we 
need to minimize 



F{w):= J2 



FT 



subjected to X]eG-B(7) ""^e = 1- Let a be the Lagrange multiplier with respect to the 
constraint J^eeEi'y) '^e- Since F is convex in We we get that F is maximized at 

maxminF(w) + a( N — 1) . (4.3) 

a w < ^ 
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So. let 

G{a) := min F{w) + > Wea 

ee£{7) 

The minimizer is obtained at 

_ 1)!!^ = ^ ^ = _ 1) W/<?|e|aV'/ 

so 

eg£;(7) 

and the minimum is obtained at 



min maxG(a)-a = maxg(g-l)i/^-iiy(^)(A)a(«-i)/«-a= (iy(^)(A))'' (4.4) 

(m,7)Gr(A) a a V / 

□ 

Proof, of Theorem 12.21 

Let us consider the sequence of reduced weighted graphs (7^, m„) (see Definition 13. 4p 
associated with a regular minimizer 7^. By Corollary I3.2l -(v) there exists a limit 7^ 
(in the sense of Hausdorff metric) of a subsequence of 7^. By (ii-iv) of the Corollary, 
\E{'yR)\ < 2N^ and is E{^r) is composed of lines. Moreover, the weights rrin '■ -E(7i?) ~^ 
M+ converges also, along a subsequence, to m : E{^r) so {m,^ji) G r(A) (see 

Definition 11.11 - (2) ) . Moreover 

lim y ml/^Je\ = V m^^lel (4.5) 

eG£;(7S) e6£;(7fl) 

By definition of the reduced graph (Definition 13. 4p and, in particular. Remark 13.11 we 
observe that X]eG-B(7g) "^^(e |e| is identical to the same expression on the non reduced 
graph 7": 

E <e\e\= E ml{!\e\. (4.6) 

eG£;{7S) eG-B(7") 
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However, on the non-reduced graphs we also have the inequahties ( 14. 1^ I4.2p 

E <e\e\<( E ^'^.^i^i') \Em\^''~'^/' = w^;\x\x-)\Em\^^-'^ 

eeE{7") \eeE(7") / 

(4.7) 

where Wg as defined in (12.21) . Here we also used the optimality of 7". 
Finally, from Theorem 12.11 

lim Wi"\A+, A-)|E(7")|(^-^)/'? = lim W^!"\x+ , X-)n'^'^-^y'' = WJX) . 

This and (I4.5M.7I) yields 

E ml/^\e\<W,{X) 

while the opposite inequality follows from the definition of Wq. □ 
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